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ABSTRACT 

Mergers of neutron stars (NS+NS) or neutron stars and stellar mass black 
holes (NS+BH) eject a small fraction of matter with a sub-relativistic velocity. 
Upon rapid decompression nuclear density medium condenses into neutron rich 
nuclei, most of them radioactive. Radioactivity provides a long term heat source 
for the expanding envelope. A brief transient has the peak luminosity in the 
supernova range, and the bulk of radiation in the UV - Optical domain. We 
present a very crude model of the phenomenon, and simple analytical formulae 
which may be used to estimate the parameters of a transient as a function of 
poorly known input parameters. 

The mergers may be detected with high redshift supernova searches as rapid 
transients, many of them far away from the parent galaxies. It is possible that 
the mysterious optical transients detected by Schmidt et al. (1998) are related 
to neutron star mergers as they typically have no visible host galaxy. 

Subject headings: gamma-rays: bursts - stars: binaries: close - stars: neutron - 
stars: supernovae 

1. Introduction 

Popular models of gamma-ray bursts (GRBs) include merging neutron stars (Paczynski 
1986, Popham et al. 1998, and references therein), and merging neutron stars and stellar 
mass black holes (Paczynski 1991, Popham et al. 1998, and references therein). However, 
the location of the recently detected GRB afterglows indicates that the bursts may be 
located in star forming regions (Paczynski 1998, Kulkarni et al. 1998a,b, Taylor et al. 1998, 
Galama et al. 1998). If this indication is confirmed with the locations of the afterglows 
detected in the future, then the NS+NS and NS+BH merger scenario will be excluded, 
as those events are expected to occur far away from the place of their origin (Tutukov & 
Yungelson 1994, Bloom et al 1998, Zwart & Yungelson 1998). 
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St ill , the mergers are certainly happening, though at the rate estimated to be several 
orders of magnitude lower than supernova rate (Narayan et al. 1991, Phinney 1991, van 
den Heuvel et al. 1996, Bloom et al 1998). It is virtually certain that a violent merger will 
eject some matter with a sub-relativistic velocity. The chemical composition of the ejecta 
must be very exotic as it is formed by a rapid decompression of nuclear density matter. It 
is not surprising that this process was suggested to be responsible for some exotic elements 
(Lattimer & Schramm 1974, 1976, Rosswog et al. 1998, and references therein). As most 
nuclides are initially very neutron rich, they will decay with various time scales. Therefore, 
we expect a phenomenon somewhat similar to a supernovae Type la, in which the decay of 
56 Ni first to 56 Co, and later 56 Fe is responsible for the observed luminosity. It is therefore 
interesting to explore the likely light curves following the NS+NS and/or NS+BH mergers. 

The neutron star mergers are expected to be among the first sources of gravitational 
radiation to be detected by the LIGO (Abramovici et al. 1992). It will be very important 
to detect the same events by other means. In the last several years a lot of effort went to 
obtain gamma-ray bursts from the mergers. However, theoretical attempts are discouraging 
(cf. Ruffert & Janka 1997, 1998, and references therein), and the observed locations of the 
burst afterglows does not favor the merger. The purpose of this paper is to point out that 
the mergers are likely to be accompanied with prominent optical transients, which should 
be detectable with the future supernova searches, and perhaps have already been detected 
by Schmidt et al. (1998). 



Modeling of a NS+NS or a NS+BH merger is very complex, and our knowledge of the 
outcome is very limited. Therefore, instead of attempting to develop a complete numerical 
model (Ruffert & Janka 1997, 1998, and references therein) we make the simplest possible 
substitute, a 'one zone model' of an expanding envelope. While most of the mass falls into 
the black hole some matter is ejected as a result of the complicated hydrodynamics of the 
merger, or the powerful neutrino burst, or the super-strong magnetic fields. For simplicity 
we assume that the expanding envelope is spherical, its mass M is constant with time, and 
its density p is uniform in space and decreases with time. The surface radius R increases at 
the fixed velocity V, i.e. we ignore dynamical effect of the pressure gradient. Therefore, the 
density throughout is given as 



2. 



Outline of the model 




(1) 
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where t is the time from the beginning of expansion. All complications of the initial 
conditions: the high temperature, the neutrino burst, and the chemical composition, are 
absorbed into several input parameters of our models: M, V, and the energy available for 
the radioactive decays. 

The temperature inside the expanding sphere varies because of several effects: adiabatic 
expansion, heat generation in radioactive decays, and radiative heat losses from the surface. 
Let us consider each of these effects. 

Adiabatic expansion. The density of expanding envelope rapidly becomes very low, 
while the injection of a large amount of heat keeps it hot, and radiation energy density 
dominates gas energy density. Therefore, we adopt 

U = 3P = aT 4 , (2) 

where U is the energy density, P is the pressure, T is the temperature, and a is the radiation 
constant. Following the first law of thermodynamics the variation of entropy S per unit 
mass is 

TdS =~ p dU+ ^Ud f ij « (^pj (fdU + 4Ut 2 dt) . (3) 

Radioactive heating. Assume that the radioactive decay of an element isotopes proceeds 
on a time scale t ra d, and releases a total amount of energy equivalent to a fraction / of the 
rest mass, so the heat generation rate per gram per second is 

fc 2 

e = - — exp (-t/t ra d) • (4) 

£rad 

If there are several decaying element isotopes, the total heat generation rate is the sum 
of that of individual element isotopes. If there are many decaying element isotopes with 
different decaying time scales, the summation can be replaced by an integration. Nuclear 
lifetimes are distributed roughly uniformly in logarithmic intervals in time, thus the total 
heat generation rate may be approximated as [] 

fc 2 

£ ~ j for ^min — t — ^max; ^min ^max- (5) 

Radiative losses. The temperature gradient is ~ T/R, the average opacity is 
k ~ k c ?s 0.2cm 2 g _1 («e is the opacity caused by electron scattering), and the radiative 



1 We are very grateful to Dr. D. N. Spergel who suggested this formula to us. 
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diffusion leads to the heat losses from the surface which are approximately given by 

where a = ac/4 is the Stephan-Boltzmann constant, T e fj is the effective temperature. 
Throughout the paper we adopt the diffusion approximation, which requires that the optical 
depth of the expanding sphere to be much larger than unity, i.e. npR 3> 1. The critical 
time t c when the expanding sphere becomes optically thin is given by upR = 1, thus 



tc 



/3kM\ 1/2 _„ , ( M V' 2 fW\- 1 f kVI 2 



\AttV 2 , 

The overall heat balance is given by 

^(e-Tf)M. ,8) 
Combining Eqs. (|1]-[8|), we obtain the equation for the variation of internal energy 

t 3 — + At 2 U « e - t 4 [7. 9 

Define 

Tl = i^L =2.80xl0 4 K' ; W M * ^ M 



3a 2 K 3 Mj \o.ooi y vo.oiM y Vk c 

(XI i t c c 



U = -^j, r=- /? = -, (10) 



then Eq. (§) can be written as 



dE7 /4 3r \ ^ _ 1 
dr I r 4/3 y r 



Z + — )^=Z5^)' (11) 



where 



t r -~ ( M \ 1/2 {W^ 1 { K\ 1 / 2 ( t r ^ 



9(r)=ae-, a = ^ « 1.13 — - H^l > (12) 

* rad l0.01M o i Ui UJ Udayi ' 1 ; 



for the case of exponential law decay; or 



9(r) = -, (13) 

T 
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for the case of power law decay. 
Define 

3/Mc 2 e o /3ttM\ 1/2 
= 4.13 x 10- erg B - [J-\ (-^f , (14 ) 

^ V - 001 / V°- O1M 0/ \«e/ 

then the luminosity of the black body radiation can be written as 

L = L q t 4 U. (15) 

The effective temperature is given by 

T cS = r l l 2 T = T^W*. (16) 
All these equations are applicable to the optically thick case, i.e. for r < 1. 



3. Solutions of the model 

There are general properties of the Eq. fliTl) which do not depend on the specific form 
of the radioactive energy generation term, g{r). At the very beginning of expansion, at 
the time of the order of a millisecond, the dimensionless time r « 1, and the internal 
energy term U is likely to be very large. The initial evolution is well approximated with the 
adiabatic expansion, and the radiative heat loss term (3rZ7/4/3), and the heat generation 
term g(r)/T 3 can be both neglected. 

The thermal evolution changes when the radioactive heat term becomes important, 
and it changes again when the radiative heat losses from the surface become important. 
However, the initial phase of adiabatic expansion makes the expanding sphere forget its 
initial thermal conditions. For this reason, in all the subsequent discussion of the analytical 
solution the initial conditions are found to be unimportant. 



3.1. Exponential law decay 

For the case of exponential law decay, g(r) is given in Eq. fll2|) . The analytic solution 
of Eq. (O) is 



-CUT 



ff Cl 4a P 6 ~ 

r 4 e 3 - 2 /8/3 + 3r 4 x 



- 6- 



x < 1 — exp 



3 ( 4 R 2 

— It ap 

8(3 V 3 H 





where Dawson's integral Y is defined by 

Y(x) = e 



''" ' e s2 ds 



-Y(-x) 



(17) 



C\ is the integral constant which is determined by the initial condition. The solution is 
practically insensitive to the initial conditions. For example, for t ~ 1 ms, T ~ 10 10 K, 
t c ~ 10 5 s, and T\ ~ 10 4 K, the solution 'forgets' the initial conditions after t = 100 s. Thus 
C\ is roughly 



3 P 



-1 +exp ( -^« 2 /^) +a \~3~ Y \ a V~3~ 



(19) 



The luminosity is 



L = L C ie - 3r2/813 + -L o a0e 



x < 1 — exp 



3 / 4 ^ 2 
— — r 

8^ V 3 H 
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Though these exact solutions are interesting for theories, they are not convenient to be 
used in practice. Thus here we also give some approximate solutions which are simple and 
convenient to be used in practice and have clear physical meanings. 

We have already defined two time scales: t rac j, on which radioactive energy source 
decays exponentially, and t c on which the optical depth of the expanding sphere is reduced 
to one. There is a third time scale td ~ npR 2 /c = f3t 2 c /t on which photons diffuse out of 
the sphere. Of particular importance is td,rad ~ o>f3t c when the photon diffusion time scale 
is equal to the radioactivity time scale, and t T fa P 1 ^ 2 t c when the photon diffusion time 
becomes comparable to the expansion time. 

If a 2 (3 ^> 1 (thus a~ l ft 1 / 2 <C a/3), the evolution of the expanding sphere can be 
divided into three regimes. The first regime is from r ~ r to r ~ (3 l l 2 . In this regime we 
may neglect the luminosity in the energy balance, i.e. we may neglect the term (3tU /4/3) 
in the Eq. (O), and we obtain 



U 
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t < L 
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L reaches maximum at r m w 1/a, and t m ~ t ra d. The second regime is from r ~ /3 1 / 2 to 
r ~ a/3. During this regime, since ar ^> 1, the radioactive term g(r) in the right hand side 
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Fig. 1. — The time variation of the bolo metric luminosity of the expanding sphere generated 
by a neutron star merger is shown for a number of models with various values of the 
logarithm of the ratio of two time scales: t c , when the sphere becomes optically thin, and 
the radioactive decay time t Ta &. The models were calculated for the fraction of rest mass 
energy released in radioactive decay / = 10~ 3 , the mass M = 10~ 2 M , and the surface 
expansion velocity V = 10 10 cm s _1 . For the adopted opacity k = 0.2 cm 2 g _1 we have 
t c = 0.975 x 10 5 s = 1.13 day, as indicated with large dots, separating the continuous 
lines (corresponding to the optically thick case) and the dashed lines (corresponding to the 
optically thin case). 
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Fig. 2.— The evolution of the model from Fig. 1 with the radioactive decay time 
^rad = t c = 0.975 x 10 5 s = 1.13 day in the effective temperature - bolometric luminosity 
diagram. The evolution begins at a very small radius, with a very high T cS and a very low 
L. The open circles correspond to the time of 0.1, 0.2, ... 1.1 days, and the vertical bar 
corresponds to the time t c when the model becomes optically thin. 
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of Eq. (IT) can be neglected, the resultant equation describes a radiation system cooling 
through expanding and radiative processes. The solution of the resultant equation is 

U - — 4 exp [-%\ , L ^ — exp ( , (22) 



«r 4 \ 8(3 

where the integral constant is determined by joining this solution with the solution before 
ft 1 ' 2 . The third regime is for r 3> aft, when the photon diffusion time scale is much 
smaller than the radioactive time scale. In such a case, the radiative energy loss is balanced 
by the energy deposited from radioactive decay, and the luminosity is simply 

L ps eM = ^-L a(3e- aT . (23) 

Since a 2 ft 3> 1, for the sub-relativistic case (ft ~ 1) this third regime happens in the 
optically thin regime, the radiation given by ( p3| ) should be non-thermal. The peak 
luminosity is L m ps Lo/a, which goes down as a increases. 

If a 2 ft <C 1 (thus a~ x ^> ft 1 ! 2 aft), the evolution of the expanding shell can be 
divided into two regimes. The first regime is from r ~ r to r ~ ft 1 / 2 . During this period, 
since «t < 1, the radioactive source g(r) = ae~ ar can be approximated by g(r) ps a. 
Inserting g(r) ~ a into Eq. (0), we get 

U « ^ (l - e" 3 ^) , L - ^Lo«/? (l - e- 3 ^ 2 ^) . (24 ) 

The second regime is for r ^> /? 1//2 , during which the diffusion time scale is much smaller 
than the radioactive time scale, the radioactive decay energy is balanced by radiative energy 
loss, and thus the luminosity is given by Eq. (p3|) . If it is extended into the optically thin 
regime, the radiation should become non-thermal for r ^> 1. The peak of the luminosity is 
roughly at the cross of the light curves given by Eqs. (|23|) and (|24"D, and the corresponding 
peak luminosity L m and time r m are 

1/2 

(25) 



4 

^aftL , 




The L m increases with a. 



From the above discussion we see that, for a 2 ft <C 1, the peak luminosity goes up as 
a increases; for a 2 ft 3> 1, the peak luminosity goes down as a increases. Thus, the peak 
luminosity reaches its maximum roughly at a 2 ft ~ 1 (or t rac j ~ ft 1 / 2 t c , which is comparable 
with t c for sub-relativistic case), where we have the biggest peak luminosity L m ~ (3 1 / 2 Lq. 



The light curves are shown in Fig. 1 for M — 10 2 M Q , V = c/3, / = 10" 
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k = 0.2 cm^" 1 (thus * c ps 0.975 x 10 5 s « 1.13 day, T x ps 2.8 x 10 4 K, L ps 4.1 x 10 44 erg/s, 
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and (3 = 1/3), and a = t c /t ra d = 10 3 , 10 2 , 10, 1, 10 _1 , and 10~ 2 respectively. The initial 
condition is taken to be T(t = 1 ms) = 2.8 x 10 10 K (but the results are insensitive to the 
initial condition). Since for (3/8(3) (r — Aa/3/3) > 1 we have L rs ^L a(3e~ aT which joins 
smoothly with L « eM, Eq. (|2"0| ) can be formally extended to optically thin region, but we 
must keep in mind that in optically thin case the radiation is non-thermal. 

The evolution of our model with the exponential decay law is shown in the log T e ^ — log L 
diagram in Fig. 2 for the case of t ra d = t c . 

The numerical results agree very well with our analytical estimates of the parameters 
important for observers: the time from the beginning of expansion to the peak luminosity, 
the peak luminosity, the temperature at the peak luminosity, and the time it takes for the 
luminosity to fall down by a factor 3 from the peak. 



3.2. Power law decay 

The exponential decay model was useful in demonstrating that the very short and 
very long time scale radioactivity is of little use for generating a large luminosity. The 
most efficient conversion of nuclear energy to the observable luminosity is provided by 
the elements with the decay time scale t ra( j comparable to t c , when the expanding sphere 
becomes optically thin. In reality, there is likely to be a large number of nuclides with a 
very broad range of decay time scales. Therefore, it is more realistic to adopt a power law 
decay model, which automatically selects the most efficient radioactive time scales. 

In the case of power law decay, g(r) is given with Eq. fllBf). The analytic solution of 
Eq. flnD is _ 

where the integration constant C4 is determined by the initial condition. Similar to the case 
of exponential law decay, in practice the solution is insensitive to the initial condition. In 
this approximation, we have C4 ~ 0. The corresponding luminosity is given by 

(27) 

For positive x, Y(x) has the maximum value at x m ~ 0.9241 and Y m ~ 0.5410. Thus, 
the time from the beginning of expansion to the peak luminosity is 

,,1/2 , ( M \ 1/2 /3V V 1/2 /'=\ 1 / 2 

K« l.5f>\ = 0.98 day {—-\ (-) (-) . (28) 
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Fig. 3. — The time variation of the bolometric luminosity of the expanding sphere generated 
by a neutron star merger is shown for models with a large mix of radioactive nuclides which 
provide a heating rate inversely proportional to time from the beginning of expansion. The 
models have three values of the fraction of rest mass released as heat: / = 10~ 3 , 10~ 4 , 10~ 5 , 
the mass M = 1CT 2 M Q , and the surface expansion velocity V = 10 10 cm s _1 . For the 
adopted opacity k = 0.2 cm 2 g -1 we have t c = 0.975 x 10 5 s = 1.13 day, as indicated with 
large dots, separating the continuous lines (corresponding to the optically thick case) and 
the dashed line (corresponding to the optically thin case). 
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log T eff [K] 

Fig. 4. — The evolution of the model from Fig. 3 with the fraction of rest mass energy 
released in radioactive decay / = 10 -3 . The evolution begins at a very small radius, with a 
very high T cff and a very low L. The open circles correspond to the time of 0.1, 0.2, ... 1.1 
days, and the vertical bar corresponds to the time t c = 0.975 x 10 5 s = 1.13 day when the 
model becomes optically thin. 



The peak luminosity is 



L m » 0.88^/% = 2.1 x 10 44 erg s" 1 ( -J-) ( — f— f ^ . (29) 
The effective temperature at the peak luminosity is 
T = „ . „ 0.79/T ^ = 2.5 x 10<K U-) (-J* ) f ^) ^ ( * ) ^ . (30) 

1 VO.ooi; v°- 01M ©/ 

From Eq. (p7[), the time when the luminosity is down a factor 3 from the peak luminosity 
is t pa 4.9/? 1 / 2 t c , which is in the optically thin epoch for the sub-relativistic case 
(V ~ 0.1c — 0.5c). For the optically thin case, the luminosity of the (non-thermal) radiation 
from the expanding shell is roughly given by L pa eM pa |L /5r _1 , which is formally the 
asymptotic form of Eq. (|2?P for r > (3 1/2 [Y(x) pa (2x)^ 1 for x 2 > 1]. Thus, Eq. (|27p can 
be formally extended to the optically thin region, where the radiation is non-thermal. The 
time from the peak luminosity to the luminosity down a factor 3 from the peak is roughly 
given by 

Fig. 3 shows light curves drawn from Eq. (|27|) with M = 1O _2 M , V = c/3, 
k = 0.2 cm 2 g _1 , and / = 10~ 3 , 10 -4 , 10~ 5 respectively. Fig. 4 shows the evolution of the 
expanding shell in the logT e ff — logL diagram for the case of / = 10~ 3 . Note that the peak 
luminosity is proportional to /, and the shape of the light curve, and all the time scales, are 
independent of / [cf. Eqs. (0-0)]. 



4. Conclusion 

Our model is so simple, that it can provide only an order of magnitude estimate of the 
peak luminosity and the time scale of a transient event that is likely to follow a violent 
merger of two neutron stars or a merger between a neutron star and a stellar mass black 
hole. The Eqs. (|27| - [31]), and Figs. 3 and 4 provide a convenient representation of our 



model, and a simple relation between the poorly known input parameters: the mass of the 
ejecta M, their velocity V, and the fraction of rest mass energy available for radioactive 
decays, and the observable parameters. Note, that for a plausible set of input parameters 
the transient reaches peak luminosity L m pa 10 44 erg s _1 , corresponding to the bolometric 
luminosity of Mbd Pa —21, i.e. in the bright supernova range. However, the duration of the 
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luminous phase is likely to be only t m ph 1 day, i.e. much shorter than a supernova. The 
duration can be extended if the ejecta have a large mass and expand slowly [cf. Eq. (ESI)]. 

While there are many improvements our model requires, the single most important is a 
quantitative estimate of the abundances and the lifetimes of the radioactive nuclides that 
form in the rapid decompression of nuclear density matter. It is possible that Rosswog et 
al. (1998) may readily provide this improvement. 

As the frequency of such events is expected to be ~ 10 3 times lower than the supernova 
rate (Narayan et al. 1991, Phinney 1991, van den Heuvel et al. 1996, Bloom et al 1998), 
they may be detected soon in the supernovae searches. The merger events are likely to be 
hotter than supernovae, which can make them easier to detect at large redshifts. The dust 
extinction which affects some supernovae is not likely to be a problem, as the mergers are 
expected far from star forming regions, many of them (perhaps most of them) outside of 
parent galaxies (Tutukov & Yungelson 1994, Bloom et al 1998, Zwart & Yungelson 1998). 

It is very intriguing that the high redshift supernova search (Schmidt at al. 1998) 
revealed mystery optical transients which typically have no host galaxy, as would be 
expected of neutron star mergers. In principle it should be fairly straightforward to test 
our suggestion by measuring future mystery events in at least four photometric bands, so 
that at least three colors would be available. In the optically thick phase our model would 
follow a locus occupied by ordinary stars in the color-color-color diagram, evolving along a 
line corresponding to a fixed redshift and changing effective temperature. As soon as the 
expanding matter becomes optically thin the colors would deviate significantly from the 
area occupied by any star with a more or less normal photosphere. 

It is a great pleasure to thank Dr. D. N. Spergel for asking us a question: 'what might 
be the observable consequences of neutron star mergers?', and for his useful comments. One 
of us (BP) acknowledges hospitality he experienced while on his sabbatical at the Institut 
d'Astrophysique, CNRS in Paris. 
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